Persistent Quantum Beats and Long-Distance Entanglement from 
Waveguide- Mediated Interactions 
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We study photon-photon correlations and entanglement generation in a one-dimensional waveg- 
uide coupled to two qubits with an arbitrary spatial separation. To treat the combination of nonlin- 
ear elements and ID continuum, we develop a novel Green function method. The vacuum-mediated 
qubit-qubit interactions cause quantum beats to appear in the second-order correlation function. 
We go beyond the Markovian regime and observe that such quantum beats persist much longer than 
the qubit life time. A high degree of long-distance entanglement can be generated, increasing the 
potential of waveguide-QED systems for scalable quantum networking. 
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One-dimensional (ID) waveguide-QED systems are 
emerging as promising candidates for quantum informa- 
tion processing [ll-Hi), motivated by tremendous exper- 
imental progress in a wide variety of systems [15l-l24j. 
Over the past few years, a single emitter strongly cou- 
pled to a ID waveguide has been studied extensively 0- 
[8l[lol[l2l4l4]. To enable greater quantum networking po- 
tential using waveguide-QED [lj, it is important to study 
systems having more than just one qubit. 

In this Letter, we study cooperative effects of two 
qubits strongly coupled to a ID waveguide, finding the 
photon-photon correlations and qubit entan glem ent be- 
yond the well-studied Markovian regime [25l-l28|. A key 
feature is the combination of these two highly nonlinear 
quantum elements with the ID continuum of states. In 
comparison to either linear elements coupled to a waveg- 
uide [29-32J or two qubits coupled to a single mode serv- 
ing as a bus [33|, both of which have been studied pre- 
viously, new physical effects appear. To study these ef- 
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FIG. 1. Schematic diagram of the waveguide system and 
single-photon transmission, (a) Two qubits (separated by L) 
interacting with the waveguide continuum. Panels (b) and (c) 
show colormaps of the single-photon transmission probability 
T and the phase shift 0, respectively, as a function of detun- 
ing 5 = ck — cjo and 2kL. Here, we consider the lossless case 
r' = 0. 



fects, we develop a numerical Green function method to 
compute the photon correlation function for an arbitrary 
interqubit separation. 

The strong quantum interference in ID, in contrast 
to the three-dimensional case (34j, m akes the vacuum- 
mediated qubit-qubit interaction [35] long-ranged. We 
find that quantum beats emerge in the photon-photon 
correlations, and persist to much longer time scales in 
the no n- Markovian regime. We show that such persistent 
quantum beats arise from quantum interference between 
emission from two subradiant states. Furthermore, we 
demonstrate that a high-degree of long-distance entan- 
glement can be generated, thus supporting waveguide- 
QED-based open quantum networks. 

Hamiltonian. — As shown in Fig. 1(a), we consider two 
qubits with transition frequencies uo\ and c^2, separation 
L = ^2—^1, and dipole couplings to a ID waveguide. The 
Hamiltonian of the system is [36| 



3=1,2 

+ ^ / dxhVj5[x-tj ) [al(x)a~ +h.c], 

j=l,2 a=R,L J 

H wg = J dx^^a t R (x)^a R (x) - a t L (x)^a L (x) , (1) 

where a^ R L (x) is the creation operator for a right- or left- 
going photon at position x and c is the group velocity of 
photons. and crj are the qubit raising and lowering 
operators, respectively. An imaginary term in the en- 
ergy level is included to model the spontaneous emission 
of the excited states at rate Y' 12 to modes other than the 
waveguide continuum [37j . The decay rate to the waveg- 
uide continuum is given by Tj = 2V^/c. Throughout the 
Letter, we assume two identical qubits: Ti = T 2 = I\ 
uoi = u 2 = ui > T. and T[ = T' 2 H. 

Single-photon phase gate. — Assuming an incident 
photon from the left (with wave vector fc), we obtain the 
single photon scattering eigenstate [38] ; the transmission 
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As shown in Fig. 1(b), there is a large window of per- 
fect transmission: T ^ 1, even when the detuning (5 = 
ck — co>o) of the single photon is within the resonance line 
width (~ r). This is in sharp contrast to the single- 
qubit case, where perfect transmission is only possible 
for far off-resonance photons Such perfect transmis- 
sion occurs when the reflections from the two qubits in- 
terfere destructively and cancel each other completely. 
Furthermore, Fig. 1(c) shows that within the resonance 
line width, there is a considerable phase shift 0. This 
feature of single-photon transmission can be used to im- 
plement a photon- atom phase gate. For example, in the 
case of 5 = — 0.5T and kL = 7r/4, the single photon 
passes through the system with unit probability and a 
7r/2 phase shift. Two successive passes will give rise to a 
photon-atom 7r-phase gate, which can be further used to 
realize a photon-photon phase gate [39| . 

Photon-photon correlation: Nonlinear effects. — To 
study the interaction effects, we develop a novel Green 
function method to calculate the full interacting scatter- 
ing eigenstates and so photon-photon correlations. We 
start with a reformulated Hamiltonian [6] 

H = Ho + V, V= £ ^d)d 3 (d)d 3 - 1), 
;=i,2 1 



H = K"3 ~ tt' 3 l2)d\d 3 + H u 



3 = 1,2 

+ E E [dxhV j 6(x-a j )[al(x)d j +h.c.], (3) 

where d}- and dj are bosonic creation and annihilation op- 
erators on the qubit sites. The qubit ground and excited 
states correspond to zero- and one-boson states, respec- 
tively. Unphysical multiple occupation is removed by in- 
cluding a large repulsive on-site interaction term U ; the 
Hamiltonians in Eqs. (pp) and (j3j) become equivalent in the 
limit U — » oo. The non- interacting scattering eigenstates 
can be obtained easily from Ho\(f>) = E\(j>). The full in- 
teracting scattering eigenstates are connected to \d>) 
through the Lippmann-Schwinger equation 111, IZol. IZl [ 
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The key step is to numerically evaluate the Green func- 
tions, from which one obtains the scattering eigenstates 
[38| . Assuming a weak continuous wave incident laser, 
we calculate the second-order correlation function g 2 (t) 
[iij for an arbitrary interqubit separation. 

Figure 2 shows g 2 (t) for both the transmitted and re- 
flected fields when the probe laser is on resonance with 



the qubit: k = k (fco = ^o/ c )- When the two qubits 
are colocated [9| (L = 0), g 2 (t) of the transmitted field 
shows strong initial bunching followed by antibunching, 
while g 2 (t) of the reflected field shows perfect antibunch- 
ing at t = 0, #2(0) = 0. This behavior is similar to that 
in the single qubit case d, Q . When the two qubits are 
spatially separated by L = 7r/2fco, we observe quantum 
beats (oscillations). Since these beats occur in g 2 (t), they 
necessarily involve the nonlinearity of the qubits and do 
not occur for, e.g., waveguide-coupled oscillators. 

As one increases the separation L, one may expect from 
the well-known 3D result that the quantum beats disap- 
pear [43]. However, in our ID system they do not: Figure 
3 shows g 2 (t) for two cases, koL = 25.57T and 100. 57r, from 
which it is clear that the beats persist to long time. The 
ID nature is key in producing strong quantum interfer- 
ence effects and so long-range qubit-qubit interactions. 

Non-Markovian regime. — To interpret these exact nu- 
merical results, we compare them with the solution under 
the well-known Markov approximation. For small separa- 
tions (koL < 7r), the system is Markovian [43[: the causal 
propagation time of photons between the two qubits can 
be neglected and so the qubits interact instantaneously. 
To understand quantum beats in this limit, we use a mas- 
ter equation for the density matrix p of the qubits in the 
Markov approximation. Integrating out the ID bosonic 
degrees of freedom yields 

dp 

dt 



^ 1 • • 



H r 
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where Fa = T -\- V while Ti2 = Tcos(ujoL/c) and Qi 2 = 
(r/2)sin(cc;o£/c) are the vacuum- mediated spontaneous 
and coherent couplings, respectively. Transforming to 
symmetric and antisymmetric states \S,A) = (\gie 2 ) ± 
\e1g2}) / V% gives a more transparent form: 
dp i r 



(3=S,A 



j3=S,A 



where cr+ A = (cr+ ± a£)/y/2, T s ,a 



(6) 



r + r ± r 12 , 

and uus,a = ^0 =•= ^12- Note that \S) and \A) are de- 
coupled from each other and have transition frequencies 
ujs,a and decay rates Ts,a which oscillate as a function 
of L. When L = 0, T s = 2T + T f and T A = T f . \S) is 
in the superradiant state, while \A) is subradiant. The 
waveguide couples only to the superradiant state and so 
the photon-photon correlation mimics that for a single- 
qubit. However, when koL = tt/2, Ts = Ta = T + T\ 
us, a = cjo ± T/2, and the waveguide couples to both \S) 
and \A). The quantum interference between the transi- 
tions \S) —> \g1g2) and I A) — » \g1g2) gives rise to quantum 
beats at frequency uos — 00 a — T, as shown in Fig.[2j 
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FIG. 2. Quantum beats in the Markovian regime. The se- 
cond-order photon-photon correlation function of both the 
transmitted (top) and reflected (bottom) fields as a function 
of t for koL = (solid line) and koL = tt/2 (dashed line). The 
incident weak coherent state is on resonance with the qubits: 
k = ko = cjo /c. (Parameters: cjq = 100r and T' = 0.1T.) 



FIG. 3. Persistent quantum beats in the non-Mar kovian 
regime. The second-order correlation function of both the 
transmitted (top) and reflected (bottom) fields is plotted as a 
function of t for koL = 25.57T (solid line) and 100. 5tt (dashed 
line). We set the incident coherent state on resonance with 
the qubits (k = k ), ljq = 100r and F f = 0.1F. 



As one increases the separation L and goes beyond the 
Markovian regime, Eq. (|5j) is not a valid description of the 
system because the causal propagation time of photons 
(or retardation effect) has to be included. Comparing the 
results in Figs. [2] and [3j we see that quantum beats are 
more visible in the non-Markovian regime in both the 
transmitted and reflected fields and persist to a much 
longer time scale, especially for the case koL = 100. 5tt. 

To better understand the persistent quantum beats, we 
extract the transition frequencies and decay rates of the 
two qubit system beyond the Markovian regime. This 
is achieved by analyzing the poles of the Green function 
(38| defined in Eq. (j4]) ; they are given by 



F(u>) 



UJ - UJQ + 



z(r + r') 



r 2 

p 2zujL/c _ q 
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(7) 



In the Markovian regime, one can safely replace uj by 
ujq in the exponent, given that ujq ^> T and L <C cT -1 . 
Eq. © then yields u± = ou - i(T + T')/2 ± iTe iuj ° L / c /2. 
The real and imaginary parts of uu± correspond to the 
transition frequencies and decay rates, which are noth- 
ing but ujs,a and — Ts,a/2 obtained using the Markov 
approximation [Eq. ©]. Beyond this Markovian regime, 
we solve Eq. © iteratively by gradually increasing L. 

Figure 4 shows that both us, A and Ts,a deviate sig- 
nificantly from their Markovian values as k^L becomes 
large [Figs. 4(c) and 4(d)]. The expanded detail plots, 
Figs. 4(a) and 4(e), show that the Markov approximation 
works well for k$L G [0,57r]. At large koL, however, both 
the symmetric and antisymmetric states become subra- 



diant \Ts,a <^ T, Fig. 4(f)]. This suppression of decay 
comes about in the following way: after the initial exci- 
tation of and emission from the first qubit, it can be reex- 
cited by the pulse reflected from the second qubit. From 
the excitation probability of the first qubit through many 
emission-reexcitation cycles, an effective qubit life time 
can be defined: it is greatly lengthened by the causal 
propagation of photons between the two qubits. Ts,a 
characterize the average long time decay quantitatively. 

The nonlinear equation Eq. $Q gives rise, of course, 
to infinitely many poles for L > 0. These poles repre- 
sent collective states of two spatially separated qubits 
with vacuum- mediated interactions. They are eigen- 
modes of the density matrix of the two qubits. The 
"two-pole" approximation of retaining only the symmet- 
ric and antisymmetric states is a good approximation be- 
cause (us, a — wq : Ts,a) are the two poles closest to the 
origin (0,0). Within the parameter range we consider, 
all other collective states are far detuned from uq and 
hence barely populated [38|. In addition, \S) and \A) 
have much smaller decay rates than all the other collec- 
tive states. Therefore, these two slowly decaying states 
dominate the long-time dynamics and quantum interfer- 
ence between their spontaneous emissions is the physical 
origin of the persistent quantum beats observed in Fig.[3l 

Qubit-qubit entanglement. — With the two-pole ap- 
proximation, we study qubit-qubit entanglement using 
the master equation Eq. ((6]) with us,a and Ts,a replaced 
by the renormalized values obtained from Eq. ([7]). We 
focus on the steady state case by including a continu- 
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FIG. 4. Renormalized transition frequencies and decay rates 
of the symmetric (S) and antisymmetric (A) states. Panels 
(a)-(c) show the transition frequencies ujs (thin solid line) 
and uja (thick solid line) obtained numerically from Eq. 
together with ujs (thin dashed line) and uja (thick dashed line) 
given by the Markov approximation. Panels (d)-(f) similarly 
show the decay rates Ys and Ya obtained both numerically 
and in the Markov approximation, (ujq = 100r and Y f = 0.ir.) 



ous weak driving laser on resonance with the first qubit: 

|28| . The entanglement is char- 
acterized by the concurrence [44]; Figure [5] shows its 
steady state value for the Rabi frequency = 0.1Y. 
For small separation [Fig. 5(a)], the concurrence agrees 
with that obtained using the Markov approximation [27j : 
C reaches its maximum when the maximally-entangled 
two-qubit subradiant state (either \S) or \A)) has a min- 
imal decay rate and is well populated [28]. Between two 
peaks, C vanishes because the symmetric and antisym- 
metric states are now barely populated and the usual 
decay rate, r + Y f > fii, holds (ia ). 

In contrast, Fig. 5(b) shows that the Markovian predic- 
tions break down: we observe enhanced entanglement for 
an arbitrary interqubit separation. Such enhancement is 
due to non- Markovian processes: both \S) and \A) be- 
come subradiant (Fig. 4) with decay rates much smaller 
than T and hence are well populated [38]. Thus, long- 
range entanglement is possible due to non-Markovian 
processes, making ID waveguide- QED systems promis- 
ing candidates for scalable quantum networking. 
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FIG. 5. Long-distance qubit-qubit entanglement. The steady 
state concurrence is plotted as a function of koL for (a) < 
koL < 57T, and (b) 95tt < koL < 1007T. The Rabi frequencies 
are Qi = 0.1T, Q2 = 0. The driving laser is on resonance with 
the qubits. (uj = 100r and Y' = OAT.) 



Discussion of loss. — Accessing the non-Markovian 
regime requires a large (effective) distance between the 
qubits and hence low loss in the waveguide. Here, we 
have included the loss of the qubit by using an effective 
Purcell factor of 10 (i.e. ~10% loss). Because waveguide 
loss has the same effect on system performance as qubit 
loss (both lead to photon leakage), we expect that the ob- 
served persistent quantum beats and long-distance entan- 
glement are robust against waveguide loss on this same 
level, namely ~ 10%. While some waveguides in current 
experimental systems are very lossy (such as plasmonic 
nanowires [15[), we can circumvent this difficulty by us- 
ing a hybrid nanofiber system as discussed in the Sup- 
plementary Material [38|. One example is an integrated 
fiber-plasmonic system [3j: the optical fiber is coupled 
to two tapered plasmonic nanowires which interact with 
local qubits (e.g. quantum dots). Another example is an 
integrated nanofiber-trapped atomic ensemble |46l HtJ: 
an optical fiber is tapered into a nanofiber in two regions 
where atomic ensembles are trapped by the evanescent 
field surrounding the nanofibers. In both of these exam- 
ples, the long waveguide connecting the two qubits is a 
high quality optical fiber in which the loss is very small 
over a length of the order of 100 wavelengths. 
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Supplementary Material for "Persistent Quantum Beats and Long-Distance 
Entanglement from Waveguide-Mediated Interactions" 

In this Supplementary Material we address the following topics: calculation of single-photon scattering eigenstates, 
our numerical Green function method, the two pole approximation, and possible low-loss systems for long-distance 
entanglement. 



Single-Photon Scattering Eigenstates 

A general single-photon scattering eigenstate of the system described by Eq. (1) in the main text reads 

\</>i) = dx (j) R {x)a) R {x) + (j) L (x)a ] L {x) + ei<r+ + e 2 cr+ |0,#i# 2 ), 



(SI) 



where |0, gig 2 ) is the zero photon state with both qubits in the ground state. The Schrodinger equation H\(j>\) — E\(j>\) 
gives 



, d 

-inc- E 

dx 



ihc- E 

dx 



b R (x) + HV 1 S(x - h) ei + HV 2 S(x - i 2 )e 2 = 0, 



(j) L {x) + HV 1 S(x - h) ei + HV 2 S(x - i 2 )e 2 = 0, 



(fiwi - iT[/2 - E) ei + Wx [Mh) + = 0, 

(hu 2 - iY' 2 /2 -E)e 2 + HV 2 [(/> R {e 2 ) + cf> L (£ 2 )} = 0. 

Assuming an incident right-going photon of wave vector k = E/c, the wavefunction takes the following form 

Akx 



(S2) 



<I>r{x) 



<j> L {x) 



2tt 

ikx 



\e{ti -x) + t 12 6(x - £i)6{£ 2 -x) + t k 6{x - £ 2 )] 



[r k 9(£ 1 -x)+r 12 6(x-l 1 )6(£ 2 -x)], 



(S3) 



where 6{x) is the step function. Setting ^^(^1,2) = [<f>R,L (£12) + <I>r,l(£i 2 )]/2 an d plugging Eq. ([S3]) into (|S2]h we 
obtain the following solution 



*12 

n 2 
tk 

Tk 



(cfc - cji + zri/2) (ck -uj 2 + iY' 2 /2 + zr 2 /2) 

(cfc - cji + ir;/2 + zri/2) (ck -uj 2 + iV 2 /2 + iv 2 /2) + f\rv^L/4' 

-*r 2 (c/c - cji + iT'j2)e 2iM * 12 
(ck - uj! + iT[/2 + zri/2) (ck — lu 2 -\- iT' 2 /2 + iT 2 /2) + T^e^/A' 
(ck - uj! + zri/2) (cfc - uo 2 + zr 2 /2) 

(cfc - cji + zri/2 + zri/2) (cfc - o; 2 + *r 2 /2 + *r 2 /2) + v x Y 2 c 2ikL i^ 

-ir 2 (cfc - cji + iT^/2 - iT 1 /2)e 2ik ^ /2 - iT^ck - u 2 + zr 2 /2 + ir 2 /2)e 2iHl /2 
(cfc - cji + iT[/2 + SV2) (cfc - cj 2 + *r 2 /2 + zr 2 /2) + r!r 2 e 2 ^ L /4 ' 

2~\ e %k£ 2 



/ 2 \ e *fe^i 



e 2 = I ?c\/ — 



(*fc — £12) • 



(S4) 



In the case of two identical qubits, tk reduces to the expression given in Eq. (2) in the main text. 

Similarly, we can solve for the single-photon scattering eigenstate for an incident right-going photon of wave vector 
k = E/c. We represent the wavefunction with an incident right-going and left-going photon by \4>i(k)) R and |0i(fc))z,, 
respectively. 
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Numerical Green Function Method 

With the Lippmann-Schwinger equation shown in Eq. (4) of the main text, we can solve for the full interacting solution. 
The non-interacting eigenstates are simply products of single-photon states. 



\<t>n(ki, * • * ,M>ai,-,a n = \<f>l (*l))ai (k 2 )) a2 ' ' " \</>l(K)) a n , ®j = R, L , j = 1-7Z, 
Ho\<f>n(ki,-' j*n))ai,-,a n = c (^l H 1" *n)|^n(*l j " " ' > *n))«i ,«» ■ 



(S5) 



For simplicity, we will focus on the two-particle solution from now on. Extending the formalism to the many-particle 
solution is straightforward. The two-particle identity in real-space can be written as 



h = Jf 



If <g> \di)(di\+I$ ®^\didj)(didjl 



i<3 



I n — ^ ^ / dx\ • • • dx n \ x\ %n) ot \ ■■■a n ('El ' ' ' X n \ i 



(S6) 



where |0) is the ground state of the two qubits (bosonic sites), \di) = d||0), \didi) = ^=^|0) and \d\d 2 ) = d\d) 2 \$). 
Inserting the above identity into Eq. (4) in the main text, we obtain 



1^2 (fcl, k 2 ))a u a 2 = \h(k 1 ,k 2 )) ai ,a 2 ^ G R (E)VI 2 \^2(kl,k 2 )) ai ,c 



\fo(ki,k 2 )) aua2 + UG R (E) MHdidilMkufo))*!, 



(S7) 



i=1.2 



Projecting Eq. ([S7|) onto (didi\ yields 



(didi |^2(^i, k 2 )) 

Oil iOCQ, 

(d 2 d 2 \ilj 2 (k u k 2 )) 

Oil ,012 



(didi 1 02 (h,k 2 )) ai>aa 
(d 2 d 2 \(/) 2 (k 1 , k 2 )) ai:OC2 



U 



Gn G12 
G 2 \ G 22 



{didi\^ 2 {k u k 2 )) 

Oil iOi<2 

{d 2 d 2 \^ 2 (ki,k 2 )) 

oil ,a>2 



(S8) 



where we introduce the short-hand notation Gij = (didi\G R (E)\djdj). Solving Eq. (|S8|) gives rise to 



{did 1 \^ 2 {k 1 ,k 2 )) Oil )Oi2 
(d 2 d 2 \^ 2 (k 1: k 2 )) aua2 



I-U 



Gn G\ 2 
G 2 \ G 22 



{didi\(j) 2 (ki,k 2 )) 

Oil ,&2 

(d 2 d 2 \(j) 2 (k u k 2 )} 

Oil ,&2 



(S9) 



Projecting Eq. ([S7|) onto a two-photon basis state (xix 2 \ and taking the U — >• 00 limit, we obtain the full interacting 
two-photon solution 



(xix 2 \M*iMaua 2 = (xix 2 \Mki,k 2 )) aua2 + U(G 1 (x 1 ,x 2 ) G 2 (x u x 2 ) ) ( ^a^l^Cifeil ife2)>«^™ 



= (xiX2\<fa(k 1 ,k 2 ))a lua2 ~G xd G d l 



(didi\(j) 2 (ki,k 2 )) 



dd I {<h<k\<h{k u k2)) au 



(S10) 



where G^^i,^) = (x\x 2 \G R {E)\didi) and 



G^d = ( Gi(xi,x 2 ) G 2 (xi,x 2 ) 

Gn G\2 
G21 G22 



Gdd 



(Sll) 



Hence, the remaining task is to calculate all the Green functions in Eq. (|S10p . This can be done using the two- 
photon non-interacting scattering eigenstates, from which we can construct a two-particle identity in momentum 
space 



oil ,o>2=R,L 

Using Eq. (|S5]h the Green functions can be evaluated as 



I 2 = ^2 dkidk 2 \(j) 2 (ki,k 2 )) 

Oil ,Oi2 

„ D T J 



(S12) 
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FIG. SI. The poles of the Green functions for (a) k L = 0, (b) k L = f , (c) k L = 25.5tt, and (d) k L = 100. 5tt. Both the 
real and imaginary parts of the poles are in units of T. We show all the poles with real part within [u>o — 4F, u)q + 4T]. The 
poles corresponding to the \S) and |j4) states are labeled as S and A, respectively. In case (d), there are four additional poles 
(Cl-CA) within the plotted range. 



dj = (didilG^E^djdj) 
= / dkidk2 

^ D T J 



1 



Gi(x u x 2 ) = {x 1 x 2 \G R {E)I' 2 \d i d i ) 



E-cki- ck 2 + i0+ 



E 



Oil ,CX2=R,L 



dk±dk2 



E — cki — ck 2 + i0+ 



(d i d i \(i)2(k 1 ,k2))a 1 ,a 2 (fa(k 1 ,k2)\d j d j ), 



(xiX2\<fa(kl,k2))a u a 2 ((h(kl, k 2 )\didi) . 



(S13) 



Doing the integrals numerically gives the full interacting two-particle solution. Again, following the same program, 
it is straightforward to extend the formalism to the three- or more photon solution with two or more qubits coupled 
to the waveguide. 

Two-Pole Approximation 

In this section, we will show the validity of the 'two-pole' approximation in the parameter regime we consider. 
Assuming two identical qubits, the poles of the Green functions in Eq. (|S10|) are given by 



i(T + T') 



a 2iujL/c 



o. 



(S14) 



Figure SI plots the poles computed numerically in four different cases. For small L, Figs. SI (a) and Sl(b) show that 
there are only two poles corresponding to \S) and \A) states within a large range of frequency. At large L, however, 
both the symmetric and antisymmetric states become subradiant [Ts,a <^ T]. This suppression of decay comes about 
in the following way: after the initial excitation of and emission from the first qubit, it can be reexcited by the pulse 
reflected from the second qubit. From the excitation probability of the first qubit through many emission-reexcitation 
cycles, an effective qubit life time can be defined: it is greatly lengthened by the causal propagation of photons between 
the two qubits. Ts,a characterize the average long time decay quantitatively. Furthermore, as L increases, there are 
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additional poles as shown in Figs. SI (c) and Sl(d), corresponding to collective states generated in non-Markovian 
processes. For L ^> cr _1 , the two-pole approximation breaks down as the additional poles of collective states become 
close enough to \S) and \A) states. 

Here, we want to analyze the case k^L = 100. 5tt [Fig. Sl(d)], where L ~ cT -1 and the two-pole approximation is 
still valid as we will show below. With a driving laser on resonance with the qubits and a Rabi frequency the 
probability to excite a state \y) {uj y , T y ) is 

p y = Z T2 VZ x2- ( S15 ) 
2+(^) 

Using this formula, we can calculate the probability of exciting the states corresponding to S (ujo + 0.32r, 0.12r), 
A (cj -0.32r, 0.12r), CI (cj + 1.08r, 0.52r) , C2 (cj + 2.0ir, 0.90r), C3 (o; + 2.98r, 1.12r) and C4 (a; +3.97r, 1.3217). 
In the limit of weak driving laser, Q —> 0, we have 

Pci = 8.6%P 5 , 
P C 2 = 2.5%P 5 , 
Pes = 1A%P S , 

Pca = 0.7%P S . (S16) 

Hence, compared to states C1-C4, S and A states are well populated and dominate the qubit-qubit interactions for 
the parameter regime considered in the main text. 

Possible Low-Loss Systems for Long-Distance Entanglement 

In this section, we discuss the issue of waveguide loss and propose several low-loss systems to overcome this difficulty. 
As discussed in the main text, waveguide loss has to be limited to the same level as qubit loss. However, some 
waveguides in current experiments, e.g. plasmonic nanowires, are too lossy to meet this criteria. We propose to use 
either a hybrid optical fiber systems or slow-light superconducting systems. In the first case, low-loss optical fibers are 
used to transmit light over a long distance. The transmission length we are considering is of order 100 wavelengths, 
thus of order 100 microns for typical quantum dots or atoms. Loss over such a distance in state of the art fiber is 
very small: taking a 4dB/km fiber, the loss will be on the order of 1 ppm. In the second case, the actual transmission 
length is very short, but due to the reduced speed of light one can still reach the non-Markovian regime. Below are 
three plausible experimental settings: (a) and (b) belong to the first case and (c) illustrates the second case. 

(a) Hybrid Fiber- Plasmonic Waveguide- Q ED System 

Figure IS2l a) shows an integrated fiber-plasmonic waveguide-QED system. The idea o f hy brid plasmonic systems was 
first proposed by Chang et al. [3]. Since then, there has been extensive experimental |49|, |50|, and theoretical [5ll-l53| 
work along this line. In the schematic, the optical fiber is coupled to two tapered plasmonic nanowires. Due to the 
subwavelength confinement [3], the plasmonic field in the nanowires couples strongly to the local qubits, e.g. quantum 
dots [IH. Coupling the nanowires to a dielectric waveguide ensures that the quantum state can be transmitted over 
long distance without being dissipated in the nanowires. 

(b) Integrated Nano fiber- Trapped Atomic Ensemble System 

In the second example, a long optical fiber is tapered into a narrow nanofiber in two regions. Then, two atomic 
ensembles are trapped by the evanescent field surrounding the nanofiber s. Strong coupling is achieved between the 
propagating photons in the nanofiber and the atomic ensembles j46[. Such a setting is a clear extension of the 
experimental systems demonstrated by several groups 46|, 47}. 



(c) Slow-light Superconducting Waveguide-QED System 

In the third example, a ID open superconducting transmission line is coupled to two superconducting qubits. It has 



been experimentally demonstrated that this system is deep in the strong coupling regime 22]. However, the typical 
length of the transmission line is on the order of the wavelength of propagating microwave photons. Hence, the 
separation of the two qubits is limited to the photon wavelength. To reach the non-Markovian regime, we can make 
the effective distance between the two qubits large by the slow-light scheme first proposed by Shen and Fan [54] . The 
idea is to couple the transmission line to an additional periodic array of unit cells made of two qubits. Flat photonic 
bands can be generated to slow down the microwave photons. While not true long-distance propagation, this could 
be an effective way to experimentally probe non-Markovian effects. 
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FIG. S2. Three possible setups for long-distance entanglement in waveguide- QED. (a) Hybrid plasmonic system. The dielectric 
waveguide (light blue) is phase-matched with the plasmonic nanowires (yellow) so that efficient plasmon transfer between them 
can be realized. The nanowires are strongly coupled to the quantum dots (blue) in the tapered regions. Note that the length 
of the dielectric waveguide between the nanowires can be very long (indicated by breaks), (b) Tapered nanofiber system. An 
optical fiber (green) is tapered into narrow nanofibers in two regions. Two atomic ensembles (red) are trapped by and strongly 
coupled to the nanofibers. (c) Slow-light superconducting system. Two superconducting qubits (blue) couple strongly to the 
slow-light superconducting waveguide (red dashed box) . 



